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Elucidating the nature of the magnetic ground state of iron-based superconductors is of paramount 
importance in unveiling the mechanism behind their high temperature superconductivity. Until re¬ 
cently, it was thought that superconductivity emerges only from an orthorhombic antiferromagnetic 
stripe phase, which can in principle be described in terms of either localized or itinerant spins. 
However, we recently reported that tetragonal symmetry is restored inside the magnetically ordered 
state of a hole-doped BaFe 2 As 2 - This observation was interpreted as indirect evidence of a new 
double-Q magnetic structure, but alternative models of orbital order could not be ruled out. Here, 
we present Mossbauer data that show unambiguously that half of the iron sites in this tetragonal 
phase are non-magnetic, establishing conclusively the existence of a novel magnetic ground state 
with a non-uniform magnetization that is inconsistent with localized spins. We show that this state 
is naturally explained as the interference between two spin-density waves, demonstrating the itin¬ 
erant character of the magnetism of these materials and the primary role played by magnetic over 
orbital degrees of freedom. 

PACS numbers: 74.70.Xa, 74.25.Ha 


INTRODUCTION 

One of the central questions to be answered in the 
iron-based superconductors is the nature of their mag¬ 
netic interactions. Because superconductivity occurs in 
proximity to a magnetic instability, it is believed that 
magnetic fluctuations play a key role in promoting su¬ 
perconducting order [HE]. In these materials, the iron 
atoms in each plane sit on a square lattice and the antifer¬ 
romagnetic state, from which superconductivity emerges, 
usually consists of stripes of iron spins aligned ferromag- 
netically along one iron-iron bond direction and antiferro- 
magnetically along the the other, with a two-fold symme¬ 
try that breaks the four-fold symmetry of the high tem¬ 
perature phase. Different theoretical approaches have 
been proposed to describe the origin of this magnetic 
two-fold (C 2 ) state, as well as the associated “nematic” 
state, and its relationship to superconductivity [3H6|. 

On the one hand, the large resistivities and enhanced 
effective masses of the iron arsenides and clralcogenides 
have been interpreted as evidence for proximity to a 
Mott transition, as seen in the similar phase diagrams 
of cuprate superconductors [7H2]- This favours an ap¬ 
proach based on localized spin models, in which the iron 
spins with fixed amplitude M, live on the sites i of 
the iron lattice and interact with each other via exchange 
interactions. This can give rise to superconductivity with 
extended s-wave symmetry, in which the order parameter 
changes sign between next-nearest neighbor sites. Some 


localized models focus not on magnetic, but on orbital de¬ 
grees of freedom, whose fluctuations in general favour a 
regular s-wave state mm • In this case, magnetic order 
is a secondary effect of the four-fold symmetry breaking 
produced by changing the relative occupation of the d xz 
and d yz iron orbitals. 

On the other hand, itinerant spin models rely on the 
metallic character of these compounds and on quasi¬ 
nesting features of their Fermi surfaces Ham]. In this 
case, instead of local spins on the lattice sites, the mag¬ 
netism is best described as a modulation of the spin po¬ 
larization of the itinerant electrons, i.e., a spin-density 
wave, S (r) = M COS (Qt), with Q — (7T, 0) or (0,7T). 
The resulting superconducting symmetry depends on de¬ 
tails of the Fermi surface, and is usually extended s- 
wave but can also be d- wave. Determining which ap¬ 
proach is valid, itinerant or localized, will therefore 
have profound consequences both for the nature of the 
emergent superconductivity and its relation to cuprate 
superconductivity [6]. 

Because both the localized and itinerant scenarios pre¬ 
dict ground states with the same space group, distin¬ 
guishing between them is a challenging task. Recently, we 
observed a new magnetic phase in hole-doped BaFe 2 As 2 
that offers a new way to resolve this issue [IT]. In pure 
BaFe 2 As 2 , the structural, orbital and magnetic tran¬ 
sitions occur simultaneously in a first-order transition. 
When sodium is doped onto the barium sites, the tran¬ 
sition temperature, T N , is reduced until magnetic or- 
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der is destroyed at x < o.3na. However, between 
0.24 < x < 0.3, four-fold (C 4 ) symmetry is restored 
inside the magnetically ordered state at T r (T r < Tjv), 
with a reorientation of the magnetic moments along the 
c-axis m • Because the magnetic Bragg peaks have the 
same reciprocal lattice indices above and below T r , this 
C 4 phase was interpreted as a double-Q magnetic struc¬ 
ture described by a coherent superposition of two spin- 
density waves, S (r) = Mi cos (Qi • r) + M 2 cos (Q 2 ■ r), 
with Qx = ( 7 T, 0) and Q 2 = (0, n). 

Although such order is naturally predicted by itiner¬ 
ant models for the iron-based superconductors [nnm- 
[ 20 ] , it was not possible to rule out the existence of an 
orbitally ordered structure that is consistent with the 
observed C 4 phase. Furthermore, the existence of mag¬ 
netic Bragg peaks at both Qi and Q 2 could be ascribed 
to domains of stripe magnetic order instead of double-Q 
order. Khalyavin et al. discuss a number of orbitally or¬ 
dered structures that are compatible with the diffraction 
m- Common to all of them is that the crystal structure 
is tetragonal, but the magnetic structure still consists of 
orthorhombic (i.e., single-Q) stripes. 

In this paper, we present unambiguous evidence 
that the magnetic state of a hole-doped iron arsenide, 
Sri_ a ,Na a; Fe 2 As 2 with x = 0.37, is a double-Q spin- 
density wave. This is the third hole-doped compound 
within the AFe 2 As 2 series (A = Ba, Sr) to show signa¬ 
tures of the C 4 phase, indicating that this is a univer¬ 
sal feature of hole-doping [3112112a- We were able to 
synthesize powders that exhibit a complete transforma¬ 
tion of the sample to the C 4 phase below T r , allowing 
us to utilize Mossbauer spectroscopy as a local probe of 
the magnetization on the 57 Fe sites. Below Tn, within 
the C 2 phase, all sites show the same Zeeman split¬ 
ting due to internal molecular fields as expected for the 
single-Q stripes. However, below the C 4 transition at 
T r < Tn, 50% of the iron sites are non-magnetic while 
the other 50% show a doubling of the magnetization, ex¬ 
actly as expected from a double-Q structure formed by 
the interference between two collinear spin-density waves, 
S (r) = Mi cos (Qi • r)+M 2 cos (Q 2 • r), with Mi = M 2 . 
This redistribution of magnetization density is not com¬ 
patible with the idea of local moments living on the iron 
sites or with any form of orbital order. Therefore, our 
results point to the primary role played by itinerant mag¬ 
netism in the phase diagram of the iron arsenides, offering 
a key insight into the nature of the electronic state from 
which superconductivity emerges. 

RESULTS 

X-ray and Neutron Diffraction 

Magnetism in the hole-doped series, Sri_ a ,Na a; Fe 2 As 2 , 
has higher transition temperatures and persists to higher 
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FIG. 1. Temperature dependent diffraction data of 

Sro. 63 Nao. 37 Fe 2 As 2 - Left panels (a-b) are from x-rays (11BM, 
A = 0.413842 A) and the right panels (c-d) are from con¬ 
stant wavelength neutrons (HB-1A, A = 2.3626 A), (a) is a 
false colour map of the data around the tetragonal (112) peak, 
and (b) is a view of the tetragonal (220) peak at temperatures 
near T r taken from the same data, (c) and (d) show the in¬ 
tensity of the (||3) and (||l) magnetic peaks, respectively, 
(e) Detailed view of the calculated intensity (green line) from 
the magnetic model fit to the 10 K POWGEN data (black 
crosses). The dotted red line shows the calculated intensities 
of a non-magnetic model. 

levels of sodium concentration than the equivalent bar¬ 
ium series [3 M- We synthesized a compound with 
the nominal composition of Sro. 63 Nao. 37 Fe 2 As 2 , below 
the critical phase boundary for magnetic order, which 
has a superconducting transition at 12 K. Rietveld re¬ 
finements of the x-ray powder diffraction spectra yielded 
a sodium concentration of x = 0.3691(5). More details of 
the sample characterization are provided in the Methods 
section and the Supplementary Information. 

The transition from tetragonal (Tt/mmm) to or¬ 
thorhombic ( Fmmm ) symmetry is evident in the pow¬ 
der x-ray diffraction as a splitting of some of the Bragg 
peaks, such as the (112) peak, shown in Figure [l^,. This 
C 2 transition, which is either weakly first order or second 
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order, occurs around Tn ~ 105(2) K, below which the or¬ 
thorhombic order parameter (*.e., the magnitude of peak 
splitting) increases rapidly with decreasing temperature. 
This behaviour is similar to many other iron-based super¬ 
conductors, but more unusually, this sample then trans¬ 
forms back to tetragonal symmetry at a strongly first 
order transition at T r ss 73 K. The first-order nature of 
this transition from C 2 to C4 symmetry can be seen in 
Figure [l]o, which shows that the two phases coexist for 
~ 10 K below T r . 

Powder neutron diffraction confirms that both the C 2 
and C 4 phases are magnetically ordered. Both the (||l) 
and (|^3) magnetic peaks are present below Tn (Figure 
[ljc and d), but the increase in intensity of the former at T r 
shows that there is a significant spin reorientation in the 
C 4 phase, as observed in Bai_ a ,Na ;r Fe 2 As 2 (0.24 < x < 
0.28) Hi]. Since the transformation back to tetragonal 
symmetry is complete in Sro. 63 Nao. 37 Fe 2 As 2 below 60 K, 
we are able to obtain a more reliable refinement of the 
magnetic structure than was possible in the earlier work. 
Figure]!]: shows that the data are fit well by a model with 
moments along c-axis, in agreement with the moment 
direction deduced by WaBer et al. PI- 

The magnetic Bragg peaks in the C 4 phase have the 
same reciprocal lattice indices as the C 2 phase, using the 
tetragonal unit cell, so one possible interpretation of the 
data is that the two phases have identical magnetic stripe 
order, only differing by the orientation of the iron spins. 
In this single-Q model, Bragg peaks from stripes parallel 
to the x and y axes in different domains would be inco¬ 
herently superposed (Fig. io- Such a model would be 
magnetically orthorhombic, so magnetoelastic coupling 
should generate an orthorhombic structural distortion as 
well, but it is plausible that it is much weaker because of 
the spin reorientation and therefore difficult to resolve. 

An alternative interpretation is that there is a single 
domain comprising a coherent superposition of magnetic 
stripes parallel to both the x and y axes, a double-Q 
model (Fig. §>)• This is the model predicted by itin¬ 
erant approaches, in which magnetic order in the C4 
phase is generated by band nesting along the x and y 
directions simultaneously, restoring tetragonal symmetry 
HU- In such a case, the residual spin-orbit coupling al¬ 
lows the parallel orientation of the resulting magnetic 
moments from each wave vector only if they are along 
the z-direction. It is well known that diffraction alone 
is unable to distinguish between multi-domain single-Q 
and single-domain multi-Q structures, since they pro¬ 
duce identical Bragg peak intensities. As discussed in 
Ref. M, it might be possible to distinguish them with 
resonant x-ray scattering, which is sensitive to the or¬ 
bital configuration of the iron d-electrons. In particular, 
space groups compatible with any possible orbital order 
would be incompatible with a double-Q model. 
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FIG. 2. Single-Q and double-Q magnetic models, (a) Single- 
Q model, in which spins are modulated with either Qi = 
(n, 0) or Q 2 = (0,7r), parallel to the a and b axes, respectively, 
in different domains, (b) Double-Q model that is formed from 
the superposition of modulations along Qi and Q 2 . 


Mossbauer Spectroscopy 

Although in reciprocal space, the single-Q and double- 
Q models look identical, they are remarkably different in 
real space. As shown in Figure (2jo, the coherent super¬ 
position of the orthogonal stripes in the double-Q model 
results in a doubling of the magnetic moments on half of 
the sites and a complete cancellation of the magnetic mo¬ 
ments on the other half. That is, half of the iron sites are 
spin-density wave nodes. In local moment systems, nodes 
represent fluctuating spins that have a high entropy, but 
in an itinerant spin-density wave, they can be a natural 
consequence of a spatially inhomogeneous magnetization. 

The best way to distinguish these two magnetic struc- 
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FIG. 3. (a-d) Mossbauer spectra (black circles) at 30, 65, 85, 
and 125 K (5, 15, 75, 95, and 105 K not shown). Fits to the 
contributions of different iron sites are shown separately: non¬ 
magnetic sites (black lines), C 4 magnetic sites (green lines), 
C 2 magnetic sites (blue lines, total spectrum (red lines), (e) 
Effective magnetic field on each magnetic site as a function of 
temperature determined from Mossbauer spectroscopy. The 
lines are guides to the eye. Error bars are smaller than the 
points. 

tures is therefore to use a local probe of the magnetiza¬ 
tion. Mossbauer spectroscopy is ideal because the Zee- 
man splitting of the nuclear levels of 57 Fe atoms is di¬ 
rectly proportional to the static magnetization density 
at the nuclear site. Earlier Mossbauer spectra on iron 


arsenides were consistent with the temperature depen¬ 
dence of the conventional antiferromagnetic stripe order, 
which is characterized by a single hyperfine field at each 
temperature [25] [25], 

We have measured Mossbauer spectra at temperatures 
between 5K and 125 K (Figure [3|. The 125 K spectrum 
shows, as expected, a single peak associated with the 
paramagnetic phase, with a small isomer shift due to 
the chemical environment that is independent of tem¬ 
perature. Spectra measured in the C 2 phase at 95, 85 
and 75 K (only 85 K is shown) are well fit with a single 
hyperfine field, characteristic of a single magnetic site, 
which grows with decreasing temperature (Figure [3^). 
However, well below the C 4 transition, at 30 K, there 
is a qualitatively different spectrum, which consists of a 
large central peak with the same isomer shift as the para¬ 
magnetic phase, indicating the presence of non-magnetic 
sites, and a sextet indicating magnetic sites with a sig¬ 
nificantly larger effective field than in the C 2 phase (by a 
factor of ~ 2). A free fit to such a two-site model shows 
that the spectral weights of each component are identical 
within the statistical uncertainty. In other words, 50% of 
the sites are magnetic and 50% are non-magnetic, exactly 
as predicted by the double-Q model. 

The 65 K spectrum, which was taken in the tempera¬ 
ture range where diffraction data indicated a co-existence 
of the C 4 and C 2 phases, shows evidence of the superpo¬ 
sition of three components, two magnetic and one non¬ 
magnetic. Although the parameters are too highly cor¬ 
related to be fit independently, the spectrum is consis¬ 
tent with a C 4 contribution, comprising an equal con¬ 
centration of large moment and non-magnetic sites, and 
a C 2 contribution from smaller moment sites. Another 
parameter, the electric field gradient, which is sensitive 
to the point-group symmetry of the surrounding ions, 
changes sign between the C 2 and C 4 phases, but is oth¬ 
erwise nearly temperature independent. The resulting 
local magnetization of the magnetic sites as a function 
of temperature (Figure^) shows a clear doubling of the 
magnetic moment within the C 4 phase compared to the 
C 2 phase, demonstrating that the C 4 magnetic structure 
involves a redistribution of magnetization density from 
the non-magnetic to the magnetic sites, an effect that is 
a clear fingerprint of an itinerant spin density wave. 

DISCUSSION AND CONCLUSION 

The Mossbauer data provide unambiguous evidence 
that the magnetism in the C 4 phase is a double-Q spin- 
density wave. This has a number of important conse¬ 
quences. As pointed out in Ref. 21] such a double- 
Q magnetic structure is incompatible with either ferro- 
orbital order involving the d xz and d yz iron orbitals or 
with more complex patterns of orbital ordering. As a 
result, it implies that the nematic phase observed in the 
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underdoped compounds is not the cause, but a conse¬ 
quence of magnetism, in agreement with the spin-nematic 
scenario [27]. Although this observation by itself cannot 
rule out the existence of orbital fluctuations, which favour 
the more conventional s-wave superconducting state, it 
does indicate the primary role played by magnetic fluc¬ 
tuations, which favour the unconventional sign-changing 
extended s-wave state. 

The most important conclusion to be drawn from this 
work is that the nature of the C 4 magnetic state is not 
consistent with a model of localized spins on the iron 
sites, in which every iron site is magnetic. At least within 
the model, widely employed as an effective model 

to study these materials 0 HSJ [29], such a non-uniform 
magnetization is not a ground state of the model. It re¬ 
mains to be seen whether modifications of this model, 
such as the inclusion of non-Heisenberg exchange inter¬ 
actions, like the biquadratic or ring exchanges, could de¬ 
scribe the non-uniform state. 

By contrast, an itinerant approach offers a natural ex¬ 
planation of this non-uniform magnetic structure as the 
interference of two nesting-related spin density waves, 
S (r) = Mi cos (Qi ■ r) + M 2 cos (Q 2 • r), with Mi and 
M 2 parallel to each other. The fact that M\ = M 2 en¬ 
sures not only the tetragonal symmetry of the system, in 
agreement with the experimental observations, but also 
implies that half of the sites are non-magnetic with their 
spin density transferred to neighbouring sites with double 
the magnetization. This is a remarkable observation that 
is only compatible with itinerant electrons. It is also con¬ 
sistent with the prediction of itinerant models that such 
a state becomes favoured over the stripe state for large 
enough doping levels [ii ini Emm ng . Furthermore, a 
secondary checkerboard charge order should accompany 
this non-uniform phase, in which the non-magnetic sites 
have locally a different charge density than the magnetic 
sites [ST]- It has been argued that fluctuations of this 
secondary charge order can enhance the extended s-wave 
transition temperature m- 

The reorientation of the magnetization along the c-axis 
follows from general group-theory arguments related to 
the space-group of a single FeAs plane with preserved 
tetragonal symmetry (32) . In the iron pnictides, spin- 
orbit coupling is not small [33], and as a consequence, 
possible spin orientations are restricted to certain crystal¬ 
lographic directions. In particular, a group-theory anal¬ 
ysis reveals three possibilities: Mi || x and M 2 || y; 
M x || y and M 2 || x; or Mi || z and M 2 || z. Because 
only the latter is compatible with the non-uniform state 
discussed here, the spins must point along the c-axis. 
This is discussed in more detail in the Supplementary 
Information. 

We note that the itinerancy of the magnetism of the 
iron pnictides does not imply that interactions are nec¬ 
essarily weak [34H36| . Indeed, even in elemental iron, a 
weak-coupling approach does not fully describe the prop¬ 


erties of the ferromagnetic state. In the iron pnictides, 
interaction effects were shown to be important to capture 
high-energy properties of the spin spectrum [37] and the 
sizable fluctuating moment observed in the paramagnetic 
state [35]. Thus, it is likely that interactions are mod¬ 
erate, and may affect distinct families of iron-based su¬ 
perconductors, such as the iron chalcogenides, differently 
00 - At least in the iron pnictides, however, our work 
demonstrates that itinerancy is an essential ingredient of 
these fascinating materials. 


METHODS 

Synthesis 

Handling of all starting materials was performed in 
an M-Braun glovebox under an inert Ar atmosphere 
(< 0.1 ppm of H 2 0 and 0 2 ). Sr (Aldrich, 99.9%) and 
Fe (Alfa Aesar, 99.99+%) were used as received. Small 
pieces of Na free of oxide coating were trimmed from 
large lumps (Aldrich, 99%). Granules of As (Alfa Ae¬ 
sar, 99.99999+%) were ground to a coarse powder prior 
to use. The precursor materials SrAs, NaAs, and Fe 2 As 
were synthesized in quartz tubes from stoichiometric re¬ 
actions of the elements at 800°C, 350°C, and 700°C re¬ 
spectively. Polycrystalline samples of Sr 0 . 67 Nao. 37 Fe 2 As 2 
were prepared from stoichiometric mixtures of SrAs, 
NaAs, and Fe 2 As, which were ground thoroughly with 
a mortar and pestle, and loaded in alumina crucibles. 
The alumina crucibles were sealed in Nb tubes under Ar, 
which were further sealed in quartz tubes under vacuum. 
The reaction mixtures were subjected to multiple heat¬ 
ing cycles between 850-950°C for durations less than 48 h 
(to minimize loss of Na by volatilization). The samples 
underwent grinding by mortar and pestle between heat¬ 
ing cycles in order to homogenize the compositions. Fol¬ 
lowing the final heating cycles, the sealed samples were 
quenched in air from the maximum temperature rather 
than allowing them to cool slowly. Initial characteriza¬ 
tion of the dark gray powders was conducted by labo¬ 
ratory powder X-ray diffraction and magnetization mea¬ 
surements. More details of the sample characterization 
are provided in the Supplementary Information. 


Powder Diffraction 

X-ray powder diffraction measurements were per¬ 
formed at Argonne National Laboratory using beamline 
11-BM at the Advanced Photon Source. Neutron powder 
diffraction measurements were performed at Oak Ridge 
National Laboratory using beamline HB-lA at the High 
Flux Isotope Reactor and the POWGEN diffractometer 
at the Spallation Neutron Source. 
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Mossbauer Spectroscopy 

Mossbauer measurements were performed in transmis¬ 
sion geometry with a sinusoidally driven 2 mCi 57 Co(Rh) 
source and a germanium detector. Silicon diode sensors 
allowed the control and stabilization of the sample tem¬ 
perature to within 0.2 K for a conventional bath cryo¬ 
stat. Powder samples having an effective area density 
of 4 mg/cm 2 of 57 Fe were placed on 99.999% pure alu¬ 
minum foil held in place by kapton tape. Calibrations 
were made using a natural a-Fe foil. The spectra were 
fit by varying the isomer shift, magnetic hyperfine field, 
and the electric quadrupole factor. The intensities of the 
magnetic sextet-split lines were constrained to a 1:2:3 
ratio according to their Clebsch-Gordon coefficients (or 
magnetic dipole matrix elements), and the Lorentzian 
linewidths for all lines of a particular iron site were con¬ 
strained to be the same. Details of the fit parameters are 
given in the Supplementary Information. 
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